#!/bin/bash
genome="hg19"
position=$1
bam_url=$2
out_wig=$3

igvtools count --bases -w 1 \
--query "${position}" ${bam_url} ${out_wig} ${genome}

